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Abstract 



We present a preliminary measurement of the central inclusive jet cross sec- 
tion using a successive combination algorithm based on relative transverse 
momenta (/cj_ ) for jet reconstruction. We analyze a 87.3 p6~^ data sample 
collected by the D0 detector at the Fermilab Tevatron pp Collider during 
1994-1995. The cross section, reported as a function of transverse momen- 
tum [pt > 60 GeV) in the central region of pseudo-rapidity (|?7| < 0.5), is 
in reasonable agreement with next-to-leading order QCD predictions. This is 
the first jet production measurement in a hadron collider using a successive 
combination type of jet algorithm. 
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I. INTRODUCTION 



Jet production in hadronic collisions is understood within the framework of Quantum 
Chromodynamics (QCD) as a hard scattering of constituent partons (quarks and gluons), 
which, having undergone the interaction, manifest themselves as showers of collimated parti- 
cles called jets. Jet algorithms associate clusters of these particles into jets in a way that the 
kinematic properties of the hard-scattered partons can be inferred and thereby compared to 
predictions from perturbative QCD (pQCD). 

Historically, only the cone algorithm has been used to reconstruct jets at hadron collid- 
ers |I|]. Although well suited to implement the experimental corrections needed in the com- 
plex environment of hadron colliders, the cone algorithms used in previous measurements 
by the hadron collider experiments at the Fermilab Tevatron present several difficulties, 
because (a) an arbitrary procedure must be implemented to split and merge overlapping 
cones, (b) an ad- hoc parameter, T^sep, is required to accommodate the differences between 
jet definitions at the parton and detector levels 0, and (c) improved theoretical predictions 
calculated at the Next-to-Next-to-Leading-Order (NNLO) in pQCD are not infrared safe, 
because they exhibit a marked sensitivity to soft radiation 0. 

Inspired by QCD, a second class of jet algorithms, which does not suffer from these 
shortcomings, has been developed by several groups These clustering or recombina- 

tion algorithms successively merge pairs of nearby vectors (partons, particles or calorimeter 
towers) in order of increasing relative transverse momentum (pt)- A single arbitrary param- 
eter, D, which characterizes approximately the size of the resulting jets, is used to determine 
when this merging stops. No splitting-merging is involved, because each vector is assigned 
to a unique jet. There is no need for introducing any ad-hoc parameters, because the same 
algorithm is applied at the theoretical and experimental level. Furthermore, by design, clus- 
tering algorithms are infrared and collinear safe to all orders of calculation. The dependence 
of the inclusive jet cross section on the choice of reconstruction algorithms or parameters is 
particularly relevant for studying the effect of hadronization and background from spectator 
partons in the event. In contrast to previous work from hadron colliders this paper 



presents the first measurement of the inclusive jet cross section using the k± algorithm to 
reconstruct jets. 



II. JET RECONSTRUCTION AND SELECTION 

The differential jet cross section was measured in bins of pt and pseudo-rapidity, t] = 
— ln[tan(6'/2)], where a right handed coordinate system is adopted with the z axis pointing 
in the proton beam direction, and 6 is the polar angle. The algorithm implemented at 



D0 1 12] is based on the clustering algorithm suggested in Ref. 0. The algorithm starts with 
a list of pre-clusters or "vectors". For each vector j, da = p^^ and for each pair of vectors, 
dij = min{p'^^,p'^ j) ARf j/D^ are defined, where D is the free parameter of the algorithm 
and ARij = ^4>ij + ^Vij is the square of the angular separation between the vectors. If 
the minimum of all da and dij is a rfjj, then the vectors i and j are merged, becoming the 
merged four-vector {Ei + Ej^ pi +Pj). If the minimum is a da, the vector i is defined as a jet. 
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This procedure is repeated until all vectors are combined into jets. Thus jets do not have 
to include all vectors in a cone of radius D, and can include vectors outside of this cone. 



The primary tool for jet detection at D0 is the uranium-liquid argon calorimeter |]T3 
which has full coverage for pseudo-rapidity |?7| < 4.1. The initial hardware trigger selected 
inelastic collisions as defined by hodoscopes located near the beam axis on both sides of the 
interaction region. The next stage required energy deposition in any A?] x A0 = 0.8 x 1.6 
region of the calorimeter, corresponding to a transverse energy (-Et) above a preset threshold. 
Selected events were digitized and sent to an array of processors. Jet candidates were 
reconstructed with a cone algorithm (with radius R = 0.7), and the event was recorded if 
any jet ii^^ exceeded a specified threshold. Jet i^^ thresholds of 30, 50, 85 and 115 GeV were 
used to accumulate integrated luminosities of 0.34, 4.46, 51.5 and 87.3 pb~^, respectively [jlO 



Jets were reconstructed offline using the k± algorithm, with D = 1.0. This value of D 
was chosen because, at next-to-leading-order (NLO), it produces a theoretical cross sec- 
tion that is essentially identical to the cone prediction for R = 0.7 |0, as used by D0 in 
its previous publications on jet production |10[. The imbalance in transverse momentum, 
"missing transverse energy" , was calculated from the vector sum of the Et values in all cells 
of the calorimeter. The vertices of the events were reconstructed using the central tracking 
system |T^. A significant portion of the data was taken at high instantaneous luminosity, 
where more than one interaction per beam crossing was probable. When an event had more 
than one reconstructed vertex, the quantity 5*^ = \Tjpj.^*\ was defined for the two vertices 
that had the largest number of tracks, and the vertex with the smaller St was retained 
for calculating all kinematic variables. To preserve the pseudo-projective nature of the D0 
calorimeter, the selected vertex was required to be within 50 cm of the center of the detec- 
tor. This requirement rejected (10.6 ± 0.1)% of the events, independent of jet transverse 
momentum. 



Isolated noisy calorimeter cells were suppressed using online and offline algorithms ||14 
Background introduced by electrons, photons, detector noise and accelerator losses that 
mimicked jets were eliminated with quality cuts. The efficiency of jet selection was approx- 
imately 99.5% and nearly independent of jet pt ■ Background events from cosmic rays or 
misvertexed events were eliminated by requiring the missing transverse energy in each event 
to be less than 70% of the pt of the leading jet. This criterion was nearly 100% efficient. 

The D0 jet momentum calibration [|12|, applied on a jet by jet basis, corrects on av- 
erage the reconstructed Pt for background from spectator partons (the "underlying event", 
determined from minimum-bias events), additional interactions, pileup from previous pp 
crossings, noise from uranium radioactivity, detector non-uniformities, and for the global 
response of the detector to hadronic jets. Unlike the cone algorithm, the k± algorithm does 



not require additional corrections for showering in the calorimeter |]T2[. For |?7| < 0.5, the 
mean total multiplicative correction factor to an observed prof 100 GeV [400 GeV] was 
1.094 ±0.015 [1.067 ±0.020]. 



III. INCLUSIVE CROSS SECTION 

The inclusive jet cross section for |?7| < 0.5 was calculated in four ranges of transverse 
momentum, using data from only one trigger in each case. The more restrictive trigger 
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was used as soon as it became fully efficient. The average differential cross section for each 
Pt bin, (i^cr/ (dpTdrj), was calculated as N/ {ArjApxe J Ldt), where Ai], Apx are the rj, pt bin 
sizes, the number of jets observed in that bin, e is the overall efficiency for jet and event 
selection, and / Ldt represents the integrated luminosity of the data sample. 

The measured cross section is distorted in prby the momentum resolution of the D0 
calorimeter. Although the resolution in jet pt is essentially Gaussian, the steepness of the 
Pt spectrum shifts the observed cross section to larger values. The fractional momentum 



resolution was determined from the imbalance in pxin two-jet events |T^. At 100 GeV 
[400 GeV] the fractional resolution was 0.061 ± 0.006 [0.039 ± 0.003]. The distortion in the 
cross section due to the resolution was corrected by assuming an ansatz function, Ap^^{l — 
'^Pr/y/s)'" , smearing this with the measured resolution, and fitting the parameters A, B and 
C so as to best describe the observed cross section. The bin-to-bin ratio of the original ansatz 
to the smeared one was used to remove the distortion due to resolution. The unsmearing 
correction reduces the observed cross section by (5.7 ± 1)% [(6.1 ± 1)%] at 100 GeV [400 
GeV]. 

The final, fully corrected, cross section for Ir^l < 0.5 is shown in Fig. ^ along with 
statistical uncertainties. Listed in Table | are the pt ranges, the point positions, the cross 
section, and uncertainties in each bin. The systematic uncertainties include contributions 
from jet and event selection, unsmearing, luminosity and the uncertainty in the momentum 
scale, which dominates at all transverse momenta. The fractional uncertainties for the 
different components are plotted in Fig. ^ as a function of the jet transverse momentum. 



IV. COMPARISON WITH THEORY 



The results are compared to the pQCD NLO prediction from JETRAD |T^, with the 
renormalization and factorization scales set to p™"^/2, where p^"^ refers to the p^of the 
leading jet in an event. The comparisons are made using parametrizations of the parton 
distributions functions (PDFs) of the CTEQ [|l^ and MRST [|l^ families. Figure |^ shows 
the ratios of {D — T)/T, where D refers to data and T to the theoretical prediction. To 
quantify the comparison in Fig. ^, the fractional systematic uncertainties are multiplied by 
the predicted cross section and a comparison is carried out. The results are shown in 
Table |T[ The agreement is reasonable {x^/dof ranges from 1.56 to 1.12, the probabilities 
from 4 to 31%), although the differences in normalization and shape, especially at low p^ , 
are quite large. The points at low p^have the highest impact on the x^- If first four 
data points are not used in the comparison, the probabilities increase to the 60 — 80% 
range. 



V. CONCLUSIONS 

In conclusion, a preliminary measurement in proton-antiproton collisions at ^/s = 1.8 
TeV of the inclusive jet cross section based on the k± algorithm has been presented. The 
quantitative test shows reasonable agreement between data and NLO pQCD predictions. 
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FIG. 1. The central (|ry| < 0.5) inclusive jet cross section obtained with the k± algorithm at 
■s/s = 1.8 TeV. Only statistical errors are included. The solid line shows a prediction from NLO 
pQCD. 
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FIG. 2. Fractional experimental uncertainties on the cross section. 
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Bin Range 


Plotted 


Cross Sec. ± Stat. 


Systematic 


(GeV) 


PT (GeV) 


(nb/GeV) 


Uncer (%) 


60 - 70 


64.6 


(8.94 ± 0.06) X 10" 


-13, +14 


70 - 80 


74.6 


(3.78 ± 0.04) X 10° 


-13, +14 


80 - 90 


84.7 


(1.77 ±0.02) X 10" 


-13, +14 


90 - 100 


94.7 


(8.86 ±0.25) X 10-1 


-13, +14 


100 - 110 


104.7 


(4.68 ±0.04) X 10-1 


-14, +14 


110 - 120 


114.7 


(2.68 ±0.03) X 10-1 


-14, +14 


120 - 130 


124.8 


(1.53 ±0.02) X 10-1 


-14, +14 


130 - 140 


134.8 


(9.19 ±0.16) X 10-2 


-14, +14 


140 - 150 


144.8 


(5.77 ±0.12) X 10-2 


-14, +14 


150 - 160 


154.8 


(3.57 ±0.03) X 10-2 


-15, +14 


1 an 1 '7r\ 


iD4.0 


/"o QQ _|_ n no\ \/ 1 n~2 
yZ.oii It KJ.yjZ) X lU 


-io, -rl'i 


170 - 180 


174.8 


(1.56 ±0.02) X 10-2 


-15, +14 


180 - 190 


184.8 


(1.05 ± 0.02) X 10-2 


-15, +14 


190 - 200 


194.8 


(7.14 ±0.13) X 10-3 


-16, +15 


ZUU — ZLK) 




on J- r> no^ in — 3 

^^4.yy it u.ucsj X iu 


-io, +io 


210 - 220 


214.8 


(3.45 ± 0.07) X 10-3 


-16, +15 


220 - 230 


224.8 


(2.43 ±0.06) X 10-3 


-16, +15 


230 - 250 


239.4 


(1.50 ±0.03) X 10-3 


-17, +16 


250 - 270 


259.4 


(7.52 ±0.23) X 10-* 


-17, +16 


270 - 290 


279.5 


(4.07 ±0.17) X 10-" 


-18, +17 


290 - 320 


303.8 


(1.93 ±0.09) X 10-" 


-18, +18 


320 - 350 


333.9 


(7.61 ±0.59) X 10-^5 


-19, +19 


350 - 410 


375.8 


(2.36 ±0.23) X 10-^ 


-20, +21 


410 - 560 


461.8 


(1.18 ±0.33) X 10-« 


-23, +27 


TABLE I. 


Single inclusive cross 


section with jets reconstructed using the k± 


algorithm in the 


central pseudo 


-rapidity region. 







PDF X^/dof Probability (%) 

MRST 26.8 1.12 31 

MRSTgt 33.1 1.38 10 

MRSTgi 28.2 1.17 25 

CTEQ3M 37.5 1.56 4 

CTEQ4M 31.2 1.30 15 

CTEQ4HJ 27.2 1.13 29 



TABLE II. comparisons (24 degrees of freedom) between JETRAD, with renormalization 
and factorization scales set to pji"-^/2, and data for various PDFs. 
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FIG. 3. Difference between data and JETRAD pQCD normalized to predictions. The outer 
error bars represent tlie total systematic uncertainty. 
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